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The central spin decoherence problem has been researched for over 50 years in the context of 
both nuclear magnetic resonance and electron spin resonance. Until recently, theoretical models 
have employed phenomenological stochastic descriptions of the bath-induced noise. During the last 
few years, cluster expansion methods have provided a microscopic, quantum theory to study the 
spectral diffusion of a central spin. These methods have proven to be very accurate and efficient 
for problems of nuclear- induced electron spin decoherence in which hyperfine interactions with the 
central electron spin are much stronger than dipolar interactions among the nuclei. We provide an 
in-depth study of central spin decoherence for a canonical scale-invariant all-dipolar spin system. 
We show how cluster methods may be adapted to treat this problem in which central and bath 
spin interactions are of comparable strength. Our extensive numerical work shows that a properly 
modified cluster theory is convergent for this problem even as simple perturbative arguments begin 
to break down. By treating clusters in the presence of energy detunings due to the long-range 
(diagonal) dipolar interactions of the surrounding environment and carefully averaging the effects 
over different spin states, we find that the nontrivial flip-flop dynamics among the spins becomes 
effectively localized by disorder in the energy splittings of the spins. This localization effect allows 
for a robust calculation of the spin echo signal in a dipolarly coupled bath of spins of the same 
kind, while considering clusters of no more than six spins. We connect these microscopic calculation 
results to the existing stochastic models. We, furthermore, present calculations for a series of related 
problems of interest for candidate solid state quantum bits including donors and quantum dots in 
silicon as well as nitrogen- vacancy centers in diamond. 

PACS numbers: 03.65.Yz; 76.30.-v; 76.60.Lz; 03.67.Lx 



I. INTRODUCTION 

The problem of decoherence of a spin interacting with 
a bath of spins (the "central spin problem" ) has its roots 
in classic works on electron and nuclear magnetic Reso- 
nance (see Ref. 1 and references cited therein). In these 
early works, the dynamics of an ensemble of spins being 
resonant with external control field (spin species A), and 
interacting with a larger ensemble of off-resonant spins 
(species B), was considered. The fluctuations of the B 
spins (due to their mutual spin-spin interactions and due 
to spin-lattice relaxation) leads to precession frequency 
fluctuations of the A spins (the spectral diffusion), which 
were then modeled as a classical stochastic process. Spin 
echo (SE) signals of A spins were calculated using differ- 
ent assumptions about the statistical properties of this 
pro cess 

The central spin decoherence problem has received re- 
newed attention due to emergence of ideas for using local- 
ized spins in solid state systems as qubits in a quantum 
computer. The currently studied systems include gate- 
defined quantum dotsP self- assembled quantum dotsf^ 
phosphorous donors in silicon,^ and nitrogen-vacancy 
(NV) centers in diamondP In all of these systems, the 
coupling of the central (qubit) spin to a bath of other 
spins is the dominant process of the loss of coherence in 



a superposition of spin up and down states (i.e., dephas- 
ing). 

A lot of attention has been recently devoted to the 
problem in which the electron spin is coupled by a con- 
tact hyperfine (hf) interaction to a bath of nuclear spins. 
For large magnetic fields only the longitudinal part of 
this interaction should be relevant (due to a large Zee- 
man splitting mismatch between the electron and nuclear 
spins suppressing their mutual flip-flops), and the deco- 
herence of the qubit should occur due to intrinsic fluctua- 
tions of the nuclear spins caused by their mutual dipolar 
coupling. Quantitative comparison between theory^^ffij] 
and experiments for spin echo in Si:P (see Refs. [T2HT5]) 
and Si:Bi (see Ref. 16) systems has shown that for natu- 
ral concentration of spinful isotope of 29 Si this is indeed 
the case. The same origin of spin echo decay was pre- 
dicted for electron spins in III-V compound base d quan- 
tum dots in the regime of large magnetic fields M. 17 l 18 l 
Recent experiments 19 in GaAs singlet-triplet qubit agree 
with these calculations for magnetic fields higher than 
~ 0.5 T. At lower fields, the electron-nuclear spin flip- 
flops cannot be completely ignored, and the SE decay is 
dominated by the contact hf interactionp2HI21 with dipo- 
lar dynamics being only a correction. 

The case of a purely dipolarly coupled system, more 
closely analogous to the original spectral diffusion prob- 
lem, was also recently brought back into focus by devel- 
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opments in spin qubit physics. One motivation is the 
fact that silicon can be isotopically enriched to reduce 
the concentration of spinful 29 Si. Below a certain concen- 
tration threshold, one can expect that the dipolar inter- 
actions between the electron spins themselves will limit 
the coherence time. In fact it was pointed out years ago 
that at large P concentrations the decay of the observed 
SE signal of the donor-bound electrons might decay due 
to dipolar interactions between these electron spins P A 
theory addressing the range of currently studied small 
concentrations of both P donors and the 29 Si nuclei was 
proposed recently! 2 ^ In that work, it was shown that (1) 
the SE decay time in Si:P is bounded by a few seconds 
due to long-range dipolar interactions between electron 
spins for realistically small donor concentrations (about 
10 13 cm -3 ) and (2) the presence of some 29 Si can actually 
increase the T 2 time considerably by suppressing donor- 
induced decoherence. The latter effect is due to the nuclei 
providing quasi-static Overhauser shifts of electron spin 
splittings, which increase the detunings between the elec- 
tron spins, and suppress the dipolar flip-flop dynamics in 
the bath! 2 - 4 ^ The predictions of Ref. 23 for SE decay times 
have been recently confirmed experimentally! 2 ^ 

Another system for which both the qubit-bath and 
the intrabath couplings are of the dipolar origin is the 
nitrogen- vacancy (NV) center in diamond PMEU j n 
case, decoherence of the qubit (which is made out of two 
levels of an electronic spin triplet) is dominated either 
by interaction with a bath of electron spins of nitrogen 
atoms (so-called PI centers), as in Refs. [27l and 129} or, in 
the case of purer samples, by interaction with a bath of 
nuclear spins of 13 C atoms, as in Refs. 26 251 andl5Tl 

In this paper we present a detailed description of a 
cluster-based theory applicable to a sparse dipolarly cou- 
pled systemf 2 ^ and we give multiple examples of applica- 
tions of the theory. In order to put this work into context, 
let us briefly review the modern microscopic approaches 
to spectral diffusion (for an attempt at pedagogical intro- 
duction to these theories see Ref. 32). A method using 
a cluster expansion of bath dynamics was developed in 
Refs. |9] and [10] and applied in the context of spin qubit 
decoherence in semiconductors. This theory produced 
results in remarkable agreement with experimental spin 
echo decay measurements using only well-known micro- 
scopic (no fitting) parameters! 3 - 3 -^ Various theories of this 
type have been applied to problems in which an elec- 
tron spin decoheres due to co ntact hf interactions with 
a dynamical nuclear spin bathP^* 17 * 33 ^ ^ In all of these 
works based upon cluster expansions of some form, the 
contributions of the bath dynamics to central spin deco- 
herence were grouped according to the number of bath 
spins participating in a nontrivial way (e.g., undergoing 
flip-flop processes). In all of these nuclear-induced spec- 
tral diffusion problems, the coupling of the central elec- 
tron spin to nuclear spins is typically much larger than 
dipolar interactions that couple the nuclear spins to each 
other. The cluster expansions are essentially perturba- 
tive expansions in the intra-bath coupling, (related to a 



diagrammatic linked cluster expansion 11 but less cum- 
bersome to compute numerically) and are therefore well 
suited to problems in which these interactions are rela- 
tively weak. Problems in which the interactions among 
bath spins are comparable to their interactions with the 
central spin (e.g., sparse, dipolar coupled electron spins) 
present a challenge for these cluster methods. In this 
article, we show that we can adapt the cluster correla- 
tion expansion (CCE) of Refs. l37l and l38l to treat these 
problems successfully. The CCE of Refs. [37] and [38] is es- 
sentially equivalent to the original cluster expansion* 9 -^ 
but greatly simplified and more convenient for consider- 
ing large cluster corrections. This method was recently 
applied to the NV center coupled to the nuclear spin bath 
where it was used to predict interesting effects related 
to the qubit back-action on the bath dynamics! 31 * 39 * 40 * 
Without relying upon the large-bath approximations or 
cumbersome corrections of the earlier cluster expansion 
theory, 10 or relying upon the clustered grouping approx- 
imation of the disjoint cluster approach,^ the CCE is 
well suited for including larger spin clusters. These larger 
clusters need to be calculated when considering dynami- 
cal decoupling of the central spi n 1 41 * 42 !, or w hen the bath 
is sparse and multi-spin correlations build within it on 
the timescale of the central spin decoherence. The lat- 
ter case applies to the problem that is the focus of this 
article. 

It may appear at first that the cluster expansion, which 
depends critically on the higher order clusters making 
systematically weaker contributions to decoherence in 
a parametrically well-behaved manner, would be com- 
pletely impractical for problems involving a sparse bath 
of environmental spins and/or a bath environment con- 
taining similar spins to the central spin. One may wonder 
that in either case (i.e., sparse bath or qubit-bath interac- 
tion being the same as the intra-bath interaction), it may 
simply be impossible to define "clusters" in any meaning- 
ful manner for a reasonable cluster expansion technique 
to work. In fact, this has inhibited the application of the 
cluster expansion technique, in spite of its great success in 
the standard spectral diffusion problem of spin decoher- 
ence in Si and GaAs, to a number of important problems 
of increasing experimental importance. In the current 
work, we establish the applicability of a cluster expan- 
sion technique for the central spin decoherence problem 
for a sparse bath with strong intra-bath couplings (com- 
parable to bath-qubit couplings). Once such a theoreti- 
cal technique is established, we can then solve a number 
of central spin quantum decoherence problems of cur- 
rent experimental relevance using it, and we apply the 
technique to solve several problems of interest in Si and 
diamond quantum computing architectures. 

The key insight, which follows from a careful imple- 
mentation of the CCE to the sparse dipolar bath and 
from extensive numerical calculations involving increas- 
ing cluster sizes, is our finding of an effect of localization 
of flip-flop dynamics of bath spins that is not obvious 
a priori. For any given group (cluster) of spins, their 
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mutual energy detunings are affected by the state of all 
the other spins, outside of the cluster. This is due to 
the diagonal (SfSj) part of the dipolar interaction. Our 
calculations show that in a sparse dipolarly coupled bath 
these interactions are introducing strong disorder in the 
energy splittings of bath spins. This disorder suppresses 
the contribution from larger cluster sizes in an appro- 
priately modified CCE. We find that the CCE can be 
adapted to converge well for a sparse dipolar bath by 
defining cluster contributions to include these externally- 
induced energy splittings and to be effectively, but effi- 
ciently, averaged over internal and external spin states. 
(We note that Ref. [17] had previously included effects of 
externally-induced energy splittings in its pair approxi- 
mation). The results and the convergence can be strongly 
dependent upon the arrangement of bath spins, but the 
convergence of results that are averaged over different 
spatial realizations of the bath is well controlled. In our 
modified formulation of the CCE, we can obtain con- 
vergent results for SE decay up to times at which the 
coherence had decayed by an order of magnitude while 
calculating clusters of at most four spins (with six spins 
clusters being shown to contribute a negligible correction 
at this timescale). 

Our main focus is thus the case of a central spin cou- 
pled to the bath spins of the same kind by dipolar in- 
teraction, e.g. an electron spin coupled to other electron 
spins. Such a situation has been extensively studied the- 
oretically by Dobrovitski et al in Refs. 1271431441 and l29l 
The most important conclusion of these papers, based 
on extensive comparison between exact numerical cal- 
culations, stochastic model, and diverse experiments for 
NV center coupled to electron spins, is that the decoher- 
ence of a single NV center (coupled to an electron spin 
bath) can be modeled very well by replacing t he bath by 
a source of classical Ornstein-Uhlenbeck noise ! 27 * 29 l 44 l For 
spin echo this means decay of the exp[— (t/TsE) 3 } form 
crossing over at long times to exp[— t/T\ ong \. An impor- 
tant distinction^ was also made between the results of 
experiments on an ensemble of qubits, and results ob- 
tained by repeated measurement of the same qubit (as 
it is done in experiments on a single NV center). This 
distinction is very important for our work here: we con- 
sider both cases (the ensemble of qubits, and a single 
qubit), since the first of them is important for current 
experiments on Si:P, while the second is relevant for cur- 
rent NV center experiments as well as considerations for 
addressable quantum bits. It is important to note that 
these studies using exact numerics are limited to very 
small bath sizes (tens of spins with current computing 
technology) and that cluster expansions do not have that 
limitation. 

The paper is organized in the following way. In Sec.[TH 
we describe the Hamiltonian of the system of interest, a 
central spin dipolarly coupled to an ensemble of randomly 
positioned spins, all of them also coupled by dipolar in- 
teractions, and define the coherence measurement pro- 
cedure (the spin echo) which our theory addresses. In 



Sec. |III[ we provide a detailed description of a variant 
of cluster expansion theory applicable to such a prob- 
lem. There we give all the details of the theory used in 
Ref. [23] to predict the 29 Si and P concentration depen- 
dence of the electron spin coherence time is Si:P system. 
In Sec. |IV[ we describe many interesting and experimen- 
tally relevant variations of the "canonical" problem de- 
fined in Sec. [Ill There we discuss the role of the intra- 
bath coupling strength relative to their interactions with 
the central spin, the possible geometrical variants of the 
problem (i.e., the case of the bath spins being localized 
in a plane some distance from the central spin), and gen- 
eralizations of spin echo experiment to sequences of mul- 
tiple pulses (dynamical decoupling). Finally, in Sec. fv| 
we present theoretical results for many example systems: 
donor-bound electrons in bulk silicon or near an interface, 
Si-based quantum dots, and NV centers in diamond. 



II. CANONICAL PROBLEM: DECOHERENCE 
IN A SPARSE ALL-DIPOLAR SPIN SYSTEM 

In this section, we describe our canonical problem of 
interest, the spin echo decoherence of a central spin in 
a sparse all-dipolar spin- 1/2 system in a strong, homo- 
geneous magnetic field environment. In this canonical 
problem, we assume that the central spin is shifted off 
of resonance from the spins of the bath, but relax this 



assumption for one of the variants in Sec. IV Using spin- 
1/2 particles is relevant for applications to electron spin 
systems; however, our methods equally capable of treat- 
ing spins of larger magnitude. 



A. System of spins 

For our canonical problem, we consider a sparse sys- 
tem of electron spins uniformly distributed at random in 
a three-dimensional continuum (or on a lattice in which 
the diluteness of the spins makes the lattice structure ir- 
relevant). The electrons are localized (bound to donors, 
for example) and dilute to the extent that they may be 
treated as point dipoles. In fact, it is not essential that 
they be electrons, but we employ conventions of nota- 
tion (such as ^-factor) that are consistent with electron 
particles. We also assume the existence of a uniform mag- 
netic field that we choose, for convention, to lie along the 
"z" direction. Our Hamiltonian, written in atomic units 
{h = 1 and l/47reo = 1), is thus 

i j>i 

where S$ are spin operators for the spin- 1/2 particles, 
Hb is the Bohr magneton, gi is the g- factor of the zth 
electron (typically, gi = 2), Bi is the externally applied 
magnetic field at each electron site, and D(r) is a tensor 
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to characterize dipolar interactions and is defined by 



A*,js(r) 



(2) 



with a,/3 = x,y,z. 5 a p is the Kronecker delta and r a is 
the a vector component of r. Our convention is to index 
the central spin as i = 0. 

For our canonical problem, we set gi = 2 and take 
the limit of a large applied magnetic field that is equal 
among all bath spins but different for the central spin. 
That is, all except the central spin are on resonance with 
each other (neglecting at this point the energy offsets due 
to dipolar interactions). If our central spin represents a 
quantum bit, it makes sense that we would be able to 
address it individually and shift its Zeeman energy to be 
off resonant from the bath spins. In taking this limit of a 
large applied magnetic field, we should disregard interac- 
tions which do not preserve the net Zeeman energy (the 
secular approximation). Thus, our effective Hamiltonian 
for our canonical problem becomes 



K 



eff ; 



i,j>0 



s: 



hz hz 



with 



3 cos 2 Si 
R ij 



(3) 



(4) 



where 6ij is the angle that the vector from spin i to spin 
j makes with the "z" unit vector (the direction of applied 
magnetic field) and is the length of this vector. By 
forcing the magnetic field of the central spin to be dif- 
ferent from the rest and taking the large field limit, we 
suppress any flip-flopping between the central spin (with 
index zero by our convention) and bath spins. This is, 
then, a standard spectral diffusion problem in which the 
polarization of the central spin is preserved but the qubit 
will dephase due to bath-induced variations of its preces- 
sional frequency (i.e., its spectral line "diffuses"). We 
later treat the case in which the central spin is resonant 



with the bath spins in Sec. IV A 

Note that the effective Hamiltonian of Eq. ([3| has 
a 1/R 3 dependence entirely. The entire Hamiltonian, 
therefore, scales with the concentration of the spins 
(Ce)- The dynamics is therefore scale invariant, with 
time that scales inversely with the Ce- Likewise, time 
scales inversely with the square of the g-factors. Our 
results, unless otherwise specified, apply to a bath of 
Ce = 10 13 /cm 3 and g = 2. However, adjusting these 
parameters only serves to rescale the time axis. 

Since we are specifying that the spins are at random 
positions in space, there are many spatial realizations of 
this problem. As a way to visualize a particular instance 
of a spatial realization of the problem, we use "celes- 
tial map" diagrams as shown in Fig. [I] These diagrams 
represent each particular "universe" from the perspec- 
tive of a central spin. Positions of the bath spins are 



projected onto a sphere (centered at the central spin) as 
cylinders whose size is proportional to the strength of 
the interaction with the central spin. The left and right 
hemispheres are split apart so we can look out in any 
direction from the central spin. We connect the repre- 
sentatives of the bath spins with rods in proportion to 
the strength of their mutual interaction as well as their 
interactions to the central spin (where these interactions 
are beyond some criteria in strength). In these way, we 
get effective "constellations" of bath spins. We chose six 
random instances labeled A through F depicted in Fig. [I] 
and will refer to these by letter throughout the text. 

We consider the limit of an infinite bath temperature in 
which the initial bath state is random without bias. We 
generate random initial bath states as product states of 
each bath spin being up or down with equal probability. 
We consider the finite-temperature variant in Sec. |IV A| 



B. Spin echo 

If one considers an ensemble of the spin systems de- 
scribed in the previous section, each with a different spa- 
tial configurations of spins and initial spin states, one 
would observe a relatively rapid dephasing occur simply 
due to the ensemble averaging. This inhomogeneous de- 
phasing time is known as T 2 *. The central spin of each 
system would experience a different shift of precessional 
frequency as a result of the magnetic field generated from 
its environment. The standard approach to remove this 
trivial effect is to apply refocusing pulses to the central 
spin. The simplest of these is the Hahn spin echo in 
which one rotates the central spin by an angle of tt about 
an axis perpendicular to the applied magnetic field mid- 
way through the evolution. A r — ^ 7r — » r sequence, for 
example, will give a refocused signal at time t = 2r. We 
report, in the study of our canonical problem, the nor- 
malized spin echo as a function of the total time, t, of 
the sequence. At t = 0, no signal is lost and we report a 
spin echo value of one. The general trend will be a decay 
of the spin echo as the pulse sequence time is increased. 
We identify the dephasing time with the time at which 
the signal reaches a value of exp (—1). 

In our study, we assume that the central spin may be 
addressed individually, and that our refocusing pulses 
are instantaneous and ideal. The methods we describe 
should be applicable to problems that relax these as- 
sumptions, and also consider other types of pulse se- 
quences, but we choose to keep the problem simple in the 
scope of this work. For ESR measurements in which all 
of the spins are essentially resonant with each other, it is 
possible to extrapolate the central spin e cho d ecay by ad- 
justing the angle of the refocusing pulse! 12 ! 25 * Otherwise, 
the decay can be dominated by the inhomogeneous de- 
cay from the environment of like spins that are all flipped 
together (known as instantaneous diffusion). By using a 
smaller angle in the refocusing pulse, the signal of the 
echo is reduced and harder to measure, but the effects of 
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FIG. 1: "Celestial map" representations of six randomly 
chosen spatial configuration instances of our canonical prob- 
lem. Each show left and right hemispheres with the applied 
magnetic field in the up (north) direction. 



instantaneous diffusion are reduced. Extrapolating to a 
refocusing pulse angle of zero yields the proper spin echo 
decay.™ 



III. METHODS 

In this section, we describe methods we use for solving 
the central spin decoherence problem with a particular 
focus on our canonical problem. We start, in Sec. |III A[ 
with a general description of the CCE method. We show, 
through examples of instances of our canonical problem, 
a need to modify the original expansion, and demonstrate 
effective techniques that overcome the arising difficulties. 



A. Cluster correlation expansion 

Explicitly evolving bath states becomes infeasible even 
for baths of moderate sizep^ with number of b ath spins 
> 20. Recently developed cluster techniques j 9 * 10 * 37 ^ 
however, can make such evaluations possible by breaking 
up the problem into smaller pieces. Our approach will 
use the CCE method. This was developed to resolve de- 
ficiencies of previously developed spin decoherence clus- 
ter expansion technique d 9 * 10 * in small bath scenerios. In 
hindsight, the formulation of Refs. |9] and [10] should be 
viewed as a large bath approximation to the CCE that 
may be convenient where applicable (in principle, the 
approximation may be systematically corrected, but if 
corrections are necessary than one is better off using the 
CCE formalism). 

The CCE has a simple and easily generalized formula- 
tion. In principle, it is always exact in the large clus- 
ter limit (apart from division by zero situations that 
may arise). In practice, the expansion converges best 
for sufficiently short simulation times but becomes nu- 
merically instable for long simulation times. Let L de- 
note the bath averaged quantity of interest. For the spin 
echo dephasing problem, we choose L = p+~ (£)/ p+~ (0), 
where Pq~(t) is the off-diagonal component of the re- 
duced density matrix for the central spin after evolving 
a t = 2r spin echo sequence. Since it is not feasible 
to solve this directly for a moderately-sized bath, let us 
define Lg, where S is any subset of bath spins, as the re- 
sult of p+~ (t) I ' p+~ '(0) as computed when we only involve 
spins outside of set S in a trivial manner (to be explained 
below). Here we are choosing to be more general than 
the original derivation of the CCE and only require that 
Ls = L when S includes the full set of bath spins and 
leave some flexibility in the way that we define Ls for 
smaller sets. We refer to a given set of bath spins, <S, as 
a cluster although there is no requirement that the con- 
stituent spins of the set necessarily be closely spaced (or 
clustered). 

At this point, the CCE formulation will approximate 
L in terms of the Ls for various S up to some maximum 
"cluster" size and will be exact when the maximum size 
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limit reaches the size of the bath. To do this, we will 
implicitly define Ls such that 

L = l[Ls, L s = n L c = L s ]J L c , (5) 

S CCS CcS 

with products over all subsets meeting the specified crite- 
ria. There will be, therefore, overlapping clusters. These 
are not disjoint sets as they are in the approach of Ref . l36l 
Explicitly, Ls is then 



L s = L s / JJ Lc 



(6) 



CCS 



as a recursive definition for any Ls- This is simply a 
tautology that serves as the definition of Ls. It becomes 
useful when we can disregard the vast majority of the Ls 
factors for L such as limiting the cluster size. We define 
the kth order of the CCE as the approximation of L with 
a maximum cluster size k: 



L 



(k) 

CCE 



n I*. 



(7) 



\\S\\<k 



Even though we are including overlapping sets of clus- 
ters, overcounting of contributions are systematically cor- 
rected as increasingly large clusters are included in the 
approximation. We know this simply by the fact that the 
CCE is exact when all clusters are included. 

To get an intuition for why the CCE takes the form 
of a product of contributing factors, consider the ide- 
alized scenario with our effective canonical Hamiltonian 
[see Eq. ([3|] in which there are two sets of non-central 
spins, A and 6, with no cross interactions (e.g., two clus- 
tered groupings that are far apart with negligible inter- 
actions to each other): 

hjeA 

+ b itj St§r-2b itj §;%, (8) 

In this scenario, with L defined as (t)/ (0), we 
may factorize L as L = Lj± x L^. Beyond the ideal sce- 
nario, this factorization is only approximate but is related 



to the perturbation theory discussed in Sec. IIIB 



We report various spin echo results calculated using 

(k) 

the CCE up to various maximum cluster sizes: L K C ^ E . 
In practice, we do not generally include all clusters of a 
given size in the calculations. We use heuristics with cut- 
off parameters to select the clusters of the most potential 
importance. We adjust the cut-off parameters until we 
are quite confident in our results. Our cluster sampling 
heuristics are described in Appendix |A{ 



B. Perturbation theory and the cluster correlation 
expansion 

(k) 

The fact that £qce ec L ua l s L in the large k limit is 
apparent from our derivation above. But how well does 




FIG. 2: Example of a connected cluster of Lemma 1 where 
edges represent the existence of bij factors of a given term of 
L c . 



this expansion converge? The key to understanding the 
convergence properties of CCE is to understand the prop- 
erties of Lc with respect to a perturbation in the interac- 
tion that couples the bath spins. Using perturbation the- 
ory, one may express Ls as an infinite power series with 
respect to the coupling constants between bath spins, de- 
noted with bij = bjj (in reference to the Hamiltonian of 
Sec. |4j but is general for any pairwise interaction Hamil- 
tonian and may be generalized for n-way interactions). 
It is also possible to expand Lc into such a power series 
by recursively expanding the denominators of Eq. 
This expansion is possible (with only positive powers of 
the coupling constants) as long as L<s, for any <S, is non- 
zero when all coupling constants are taken to be zero; in 
the case of L = p+~(t)/p+~(0) and using our effective 
Hamiltonian in Eq. Ls = 1 + 0(bijt) and, by the 
recursive definition of Eq. ([H]), Lc = 1 + 0(bijt). 

As defined and noted in Ref. |37j L c = 1 + 0{b^ j t k ) 
with k = ||C||. In our generalization, such a result is con- 
ditional. The following definition, lemma, and theorem 
specify the conditional perturbative properties of Lc. 

Definition 1 

We say that Ls is factorable under disconnected in- 
teractions if, for any S, given X G S, y G S such 
that X [J y = S , X C\ y = $ (disjoint), and that all 
biexjey = °; then L s = L x L y . 

Lemma 1 // Ls is factorable under disconnected inter- 
actions, all non-constant (bij dependent) terms of Lc 
in a power expansion with respect to bij, must contain 
factors of bij that, when viewed as graph edges between 
nodes i and j, connect all spin "nodes" in C fully. In 
other words, each non-constant term will involve all spins 
in C and may not be factored into parts that involve dis- 
joint, non-empty sets of spins (Fig.^ty. 



Theorem 1 If Ls is factorable under disconnected in- 
teractions, Lc = const. + 0(bi 
(the size ofC). 



1 t k ~ 1 ) where k = ||C|| 



The above Theorem follows directly from the Lemma 
considering the simple fact in graph theory that a set 
of k nodes cannot form a connected graph with less than 
k — 1 edges; the Lemma is proven in Appendix [Bl The 
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lowest non-constant order will contain an additional fac- 
tor of bij under certain circumstances, so that Lq = 



const. 



0(bijt k ), but this is not an important consider- 
ation for this discussion. 

The important point is that larger clusters will intro- 
duce higher order corrections with respect to bij. This is 
the essence of the CCE^and previou s 9 * 10 * cluster expan- 
sions. However, bijt <C 1 is not a strict requirement for 
convergence; a localization effect in disordered systems 
has been demonstrated to extend convergence into the 
bijt > 1 regime^. This in fact happens in our canoni- 
cal problem, albeit a few modifications of the CCE (ex- 
plained in detail below) are necessary to capture this ef- 
fect. In any case, the best practice is to test convergence 
by computing an extra order in the expansion, providing 
an estimate of the error. 

The CCE works extremely well out into the tail of the 
spin echo decay in those central spin decoherenc e pro b- 
lems, such as nuclear-induced spectral diffusion,™^ in 
which the central spin is strongly coupled to many bath 
spins relative to the coupling strength among bath spins. 
In this case, the decay timescale is typically short com- 
pared with interaction strength among the bath spins 
(i.e., the bath is slowly evolving). Our canonical prob- 
lem pushes us into the challenging regime of this pertur- 
bation because the central spin has interactions to the 
bath spins that are comparable to the strength of the 
interactions among the bath spins. This presents chal- 
lenges to the cluster method, but we shall demonstrate 
effective techniques to address these challenges. 



C. Treatment of "External" Spins 



The CCE that we described in Sec. lIII Al has some flex- 
ibility. We have stated that L$ is to be defined in a man- 
ner that involves external spins trivially. The most simple 
way to define L$ would be to ignore all the external spins 
entirely. That is, L$ is the result of (t) / p\~ (0) when 
only including those bath spins into the problem that are 
contained in set S. We show the CCE results from apply- 
ing this definition for the spatial configuration instance 
A (Fig. B in Fig. R At short times, the L^ce (up to 3- 
clusters) exhibits a small correction to L K C ^ E (2-clusters). 
The CCE, in this form, does not provide a robust solution 
to our canonical problem except to study the initial part 
of the decay. Before significant decay occurs, L^ E blows 
up and becomes numerically unstable. We note that the 
primary contribution from 3-clusters near the onset of 
this numerical instability comes from the suppression of 
flip-flopping dynamics due to the magnetic field gradi- 
ents generated by external spins as depicted in the upper 
cartoon of Fig. [3] That is, the 2-cluster contributions are 
overestimated when we completely ignore all of the ex- 
ternal interactions, and this must be compensated at the 
3-cluster level. When we ignore all other bath spins, any 
pair of spins is completely resonant and can freely flip- 




10 

Time (s) 



10 

Time (s) 



FIG. 3: (Top) Depiction of a 2-cluster whose flip-flops are 
suppressed by the magnetic field gradient generated by a third 
spin. (Bottom) Spin echo results from the CCE applied to 
spatial configuration instance A (see Fig. [TJ when we don't 
define a cluster contributions to be "externally aware." Each 
"spaghetti" strand is the result for a different random initial 
spin state (a product state of up or down for each spin) with 
the mean of L^ CE [L^ E ] as encircled +'s [triangles]. The 
3-cluster contributions must compensate for the lack of "ex- 
ternal awareness" which leads to numerical instability at later 
times. 



flop. Considering the presence of these other bath spins, 
they are generally off-resonant to some degree. That is, 
the Ising-like part of the long-range dipolar interactions 
plays an important role regarding whether or not a given 
pair of spins, for example, are resonant with each other 
for flip-flopping. Capturing this effect for small clusters 
is crucial for obtaining our computationally feasible and 
convergent theory of spectral diffusion in the canonical 
problem. 

Ideally, one would want to define things in such a way 
that Lemma 1 would be applicable with respect to just 
the flip-flop interactions and that the SfSj interactions 
would come for free. We have not found such an effi- 
cient solution that achieves this, but a step in the right 
direction is to include these Ising-like interactions with 
spins outside of S for a given Ls- That is Ls excludes 
only the flip-flop interactions involving spins external to 
S. For each bath state |J), let L J = (J+ (t)\J~ (t)) 
where J ± (t) = U(t)\ J, ±), using ± for the up/down cen- 
tral spin states. Now define L J S = (J$ (t)\Jg (t)) where 
\Jg(t)) = Us(t)\J, ±); here, Us gives the evolution which 
disregards all except the Ising-like SfS? interactions with 
bath spins external to S. We thus have the CCE defined 
for each J state independently. For a given density ma- 
trix p = J: jP j\J)(J\, L = EjPjL j = (L J s )j. 

We show the results of these CCE calculations, with 
"external" spin awareness as described above, in Fig. H] 
for each of the spatial configurations depicted in Fig. [l| 
Comparing with Fig. |3j considerable improvement is ap- 
parent. However, these results, which still show strongly 
unstable behavior at longer times can be further im- 
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proved. 



D. Interlaced spin state averaging 

As we have just discussed, disorder in a system may 
improve the CCE convergence due to a localization. It 
can also be problematic and create numerical instability, 
however. The problem arises when some Lq attains, by 
chance of the disorder, a very small value for a time that 
is short compared with the overall decoherence time. In 
the calculation of L$ for some supercluster S D C, Lq 
may be a factor in its denominator so that L$ attains a 
very large value, disproportionate to its order, i.e., ||<S||, 
in the CCE. We note this effect in the numerical insta- 
bility observed in Fig. [4] for later spin echo times. 

We find that this ill effect from disorder may be mit- 
igated very effectively by defining the CCE to use "in- 
terlaced" spin state averaging. Rather than computing 
a separate L for each given bath state J, we can self- 
consistently average over the bath states internal to the 
computation of L$- In terms of the definitions of the 
previous sub-section, we may define Ls = (Lg)j, as an 
average over bath states J, and use the standard CCE 
equations [see Eqs. ([5|, ([6]), and ([7])]. In the large cluster 
limit, this will approach the correct solution. Defined in 
this way, Ls is not factorable under disconnected inter- 
actions (this is because (S z ) 2 ^ ((S z ) 2 )), so the lemma 
and theorem of Sec. Ill B| does not carry through. How- 
ever, we find the convergence behavior to be significantly 
improved (see Fig. [6| by using this sort of strategy. 

It is only feasible to approximate Ls with this spin 
state averaged definition because explicit averaging over 
the state of all external spins has exponential complex- 
ity. In practice, we instead use the following formulation 
which approximates interlaced spin averaging in a man- 
ner that corrects itself as more clusters are included in 
the approximation. Let T be a set of clusters (e.g., up to 
a certain size) that we include to approximate the solu- 
tion. Let J be some bath spin state, as a product state 
of up or down for each spin, that will serve as a template. 
We now define 



2-cluster only 



up to 3 -clusters 



(9) 
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where /C(J, C,T) is the set of all spin states that may 
differ from J only for spins in superclusters of C that are 
contained in T. That is, 

/c(j,c,r) = 3C e r,C d c,p(|j), c c'}, 

(10) 

where V{\ J), | J')) is the set of spins whose state differs 
between \J) and \J'): 



FIG. 4: ^cce (l e ft) and £qq E spin echo results correspond- 
ing to spatial configuration instances A-F (Fig. [TJ of our 
canonical problem using the "externally aware" treatment of 
clusters. Each "spaghetti" strand is the result for a different 
random initial spin state (a product state of up or down for 
each spin) with the mean of L^ CE [£qq E ] as encircled +'s 
[triangles] . 



V 
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2-cluster only 



up to 3-clusters 



up to 4-clusters 



FIG. 5: Illustration of /C(J,C,T) [see Eq. (To)]. Bath spins 
are denoted by the circles; up/down arrows denote the spin 
state "template" J. Clusters in F (the included clusters) are 
denoted by ovals that encompass multiple bath spins, C be- 
ing the one with the solid, thick outline. Superclusters have 
dashed outlines and all other clusters have dotted outlines. 
/C(J, C,r) are all spins states, including J, that only differ 
from J in spins within superclusters of C, shown with the 
thicker outline and t^4- or inscriptions (the up/down 

arrow on the left, say, denotes the state of the J template, 
but all combinations of these are included in JC). 



An example for /C(J, C,T) is illustrated in Fig. |5j Then 
we define 



(12) 



C'dC 



where L J C solves the Lc problem for the given spin state 
J. Importantly, this yields the exact spin state average 
solution for C E in the limit that T includes all clusters 
(J becomes irrelevant). Furthermore, it may be com- 
puted relatively efficiently if we are sufficiently selective 
in choosing the clusters of T (see Appendix |A|). With 
proper bookkeeping, each Hamiltonian (for a given clus- 
ter and external spin state) need only be diagonalized 
once, and each L J C need only be computed once and raised 
to the proper power to be multiplied into the solution. 

We present results using this revised CCE for our six 
respective spatial configurations (Fig. [I]) in Fig. [6j The 
numerical instabilities we observe in Fig. [4] have been re- 
moved, and the convergence appears to be well behaved 
as we go from L^ E on up to L^ CE . We still see some 
unphysical results (larger than unity) for some spin state 
templates, but not as widely ranging and erratic as be- 
fore without interlaced spin averaging. In principle, the 
CCE results should be exact when all clusters are in- 
cluded and any non-physical results would go away. This 
trend toward the physically valid range is oberved in go- 
ing from Lq£ E to Lqq E (Lq^ e always gives results in 
the physical range). We also note, as an indication of 
convergence, that a split between averaged 4-cluster and 
3-cluster results in each case (particularly visible in spa- 
tial configurations A, C, and F) occurs later in the time 
parameter than the split between averaged 3-cluster and 
2-cluster results. The convergence is not always great for 
this challenging problem, but a good fraction of the de- 
cay appears to be captured well. As an aid to intuition 
about what is going on for the different cluster sizes, we 
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FIG. 6: L 

echo results corresponding to spatial configuration instances 
A-F (see Fig. [I]) of our canonical problem using "external 
awareness" and "interlaced spin averaging" in our implemen- 
tation of the CCE as expressed in Eq. Each "spaghetti" 
strand is the result for a different random spin state template 

J [see Eq. j^]. The mean of L^ CE , £(?CE' an< ^ ^CCE are re P" 
resented as encircled +'s, triangles, and squares respectively. 
The brown (color online) dashed curves of the form exp (— t 3 ) 
on the right panels are presented for comparison. 
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FIG. 7: Depiction of important 3-cluster (top) and 4-cluster 
(bottom) contributions in our adapted CCE with "external 
awareness" and "interlaced spin averaging." At the 2-cluster 
level, all flip-flopping pairs are treated independently. The 
3-cluster level must compensate for overlapping flip-flopping 
pairs that cannot be independent because they share a bath 
spin. The 4-cluster level must compensate for the fact that 
two separate flip-flopping pairs in close proximity influence 
each other magnetically (through the Ising-like SfSj iterac- 
tions) as they flip-flop. This can enhance spectral diffusion, 
for example, if the two pairs are off resonance independently, 
but on resonance (conserve energy) together as a simultane- 
ous process. 



depict important 3-cluster and 4-cluster contributions in 
Fig.0 

We include decay curves of the form exp (— t 3 ) on the 
right panels of Fig. [6] for comparison with the calcula- 
tions. Such behavior is expected in the initial decay 
for Ornstein-Uhlenbeck noise. With exception to spatial 
configuration A, the exp(— t 3 ) form fits the calculated 
results very well for the first 25% to 50% of the decay, 
confirming the results from Refs. [27j [44j and [29] using 
a calculation starting from a microscopic model of the 
dipolarly coupled bath. 



E. Ensemble of spatial configurations 



2-cluster only up to 4-clusters 




time (s) time (s) 

FIG. 8: ^cce (l e ft) and £qq E (right) spin echo results of 
our canonical problem using "external awareness" and "inter- 
laced spin averaging" in our implementation of the CCE as 
expressed in Eq. We use the scaling parameter f: time 
is scaled inversely with £ and q 2 Ce = ( x 4 x 10 13 /cm 3 (i.e., 
C = 1 is the result for g = 2 and C E = 10 13 /cm 3 ). Each 
"spaghetti" strand is the result for a different random spatial 
configuration averaged over a large number of spin state tem- 
plates. The mean of £(?CE an< ^ ^CCE are re P resen ted as en- 
circled + 's and squares respectively. The median of L^ CE and 

^CCE are re P resen ted as encircled x's, and diamonds respec- 
tively. The bottom figures show the same spin echo results 
on a logarithmic scale for the decay. 

(two, four, and six), we get slight corrections that push 
down the tails of the spin echo decay curves. Let us stress 
again that the fact that only relatively small clusters of 
flip-flopping spins are enough to describe the SE decay is 
due to our use of a modified version of CCE. Averaging 
of these cluster contributions over the states of external 
spins, while keeping the diagonal dipolar interactions be- 
tween the spins from the cluster and from its outside, 
allows us to capture the effect of localization of spin dy- 
namics over the duration of appreciable coherence. 



We show results for a large ensemble of different spa- 
tial configurations in Fig. [8] We indicate the median and 
mean results for the ensemble (well converged to approx- 
imate the infinite ensemble). Both the median and mean 
are computed separately for the statistical results at each 
point in time; that is, the median values at different time 
points do not necessarily come from the same spatial con- 
figuration result. The median and mean differ drastically 
at short times because the mean is dominated by rare 
cases where the spin echo dips down early (for example, 
when a pair of bath spins happen to be particularly close 
to the central spin). Our CCE expansion, with "exter- 
nal awareness" and "interlaced spin averaging" exhibits 
very good convergence as demonstrated most clearly in 
Fig. [9j as we increase the cluster size in even numbers 



F. Convergence 

We demonstrate convergence of the ensemble average 
spin echo for our canonical problem in Fig.[9j To be confi- 
dent in our results, we must ensure sufficient cluster sam- 
pling (see Appendix [A]) and the inclusion of sufficiently 
large clusters. The cluster sampling described in Ap- 
pendix [A] relies upon radial cutoffs and cluster counting 
cutoffs. Our confidence in our results increases as we in- 
crease either of the cutoffs. We employ convenient meth- 
ods that automatically increase the cutoffs as needed by 
comparing results from differing cutoff values; for exam- 
ple, if there is a significant difference in results that in- 
clude different numbers of clusters, we increase the clus- 



ter count cutoff. 

It is best to evaluate and analyze the effects of larger 
clusters in terms of their correction to solutions that ex- 
clude them: Lqq E — Lqq-^K For the most rapid con- 
vergence, this correction (difference) should be averaged 
over the ensemble of spatial realizations of the bath 

(k) 

rather than the individual i^CE resm ^s. Such correc- 
tions, for successive cluster sizes, are presented in the 
lower panel of Fig. [9] The corrections are only significant 
toward the tail of the decay with an onset time that in- 
creases with increasing cluster size. At later times, larger 
cluster corrections will often become numerically unsta- 
ble, but the ealier parts of the decay are most important 
for quantum information considerations. 

We note that odd cluster sizes tend to increase the spin 
echo decay time (enhanced coherence) and even cluster 
sizes tend to decrease the spin echo decay time (further- 
ing decoherence) . Furthermore, 4-cluster corrections are 
the same order of magnitude as 3-cluster corrections and 
6-cluster corrections are the same order of magnitude as 
5-cluster corrections. The reason for this even/odd trend 
relates the spin up / down symmetry when regarding non- 
polarized baths and was noted in Ref. [10] for the case of 
nuclear-induced spectral diffusion. A similar argument 
applies here. When averaging over a bath described by 
density operator proportional to unity (as it is the case 
at high temperatures and low magnetic fields), the oper- 
ation of reversing the sign of the Hamiltonian leaves the 
expression for spin echo signal unchanged [see Eqs. (60) 
and (61) in Ref. [TO] . Therefore, in the non-polarized av- 
erage, only even perturbation order terms survive. Also 
note that contributing processes must perform a full cy- 
cle in state space, looping back to some initial state; that 
is because the bath states are traced out in obtaining the 
reduced density matrix of L. Thus, 3-clusters are fourth 
order in such a perturbation expansion and 5-clusters are 
sixth order. For this reason, also, 3-clusters are not dom- 
inated by the process of three flip-flops since this is can- 
celed by the symmetry of the unpolarized bath [this was 
also noticed in Ref. [TTJ see the discussion there above 
Eq. (12)]. Instead they are dominated by double flip- 
flops of two pairs sharing one spin in common as in the 
top panel of Fig. [7[ This provides a correction to the 
2-cluster level which overcounts these as two indepen- 
dent flip-flopping pairs; at the 3-cluster level, we see that 
the shared spin cannot simultaneously flip-flop with two 
partners. Presumably, a similar situation occurs at the 
5-cluster level, though the picture is substantially more 
complicated. As a side note, a nice aspect of the CCE 
over the linked cluster expansion^ in practice is that we 
do not have to decompose this complicated picture to 
obtain results. 



IV. VARIATIONS OF THE PROBLEM 

Our canonical problem was chosen for its simplicity as 
a problem to study in depth without complicating details 
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FIG. 9: Ensemble- averaged spin echo results and correc- 
tions for our canonical problem within various approxima- 
tions using our adaptation of the CCE. The upper panel gives 
^CCE f° r var i° us cutoffs in our cluster sampling heuristics 
(Appendix [Aj. The £ scaling parameter is used as in Fig. [8] 
Relative to corresponding solid line curves, dotted lines in- 
clude twice as many clusters (Nk), dashed lines increase the 
radial cutoff (Rc) by 30%, and dot-dashed lines do both. The 
blue curves (color online) in the top panel demonstrate cut- 
offs that are insufficient at later time scales (N2 — 200 and 
Rc = 600 nm for the solid blue line). The black curves use 
good cutoff values (N k ~ 30,000 and R c ~ 1000 nm); for 
this reason, it is difficult to distinguish corresponding curves 
of differing line patterns. The lower panel shows successive 
corrections as we increase the maximum cluster size. 

to serve as distractions. In the study of real- world prob- 
lems, many variations of our canonical problem arise. 
In this section, we consider some natural variations of 
this problem and discuss resulting trends in generality. 
The canonical problem will provide a convenient refer- 
ence point as we analyze these trends. We will reference 
each of these variations in Sec. |V| where we discuss spe- 
cific problems of interest in the realm of (solid state spin) 
quantum information. 



A. Coupling strengths, polarization, and 
resonances 

In this section, we consider variations of the canonical 
problem with regard to coupling strengths, resonances 
among pairs of spins, and polarization of the bath spins. 

We first consider the effects of adjusting the coupling 
strengths by simply adjusting the g-factor of the central 
spin. Only the relative strength of intra-bath interactions 
versus central spin-bath spin interactions are of qualita- 
tive importance (the absolute strengths only affect the 
scale of the time parameter). By adjusting go and hold- 
ing ^>o constant, we explore the transition from a cen- 
tral spin that is weakly coupled to the bath (go small) 
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FIG. 10: Comparison of Hahn echo results for different go 
values (the central spin g-factor) in the large magnetic field 
limit and off-resonant central spin treatment (e.g., the canon- 
ical problem otherwise). The £ scaling parameter is used as 
in Fig. HI with g^oCs = C x 4 x 10 13 /cm 3 . The upper panels 
display £qq E (+'s), ^cce (triangles) and ^cce ( sc t u ares) 
mean value results. The lower panels display ^^CE median 
(diamonds) as well as mean (squares) value results on a log- 
arithmic scale. 



to one that is strongly coupled to the bath (go large). In 
relative terms, this is also an exploration of the transi- 
tion from strong (g^o large) to weak (#i>o small) cou- 
pling within the bath respectively, although this requires 
a translation of our timescales to reflect holding go con- 
stant instead of holding g^o constant. We again treat 
the case in the large magnetic field limit where the central 
spin is off-resonant with the bath spins (either by a local 
magnetic field offset or the g-factor difference), as char- 
acterized by Eq. (J3|. We compare three difference cases 
in Fig. [ToJ go = ^>o/10 (weak coupling to the bath), 
9o = 9i>o (comparable interactions), and go = 10^>o 
(strong coupling to the bath). The time scaling factor 
is dependent upon the value of gi>o- Increasing go in- 
creases the coupling of the central spin to the bath with 
causes a decrease in decoherence time. It also makes 
sense that larger clusters become more important in the 
regime in which the coupling within the bath is relatively 
strong. In our go = ^>o/10 case, the cluster expansion 
is only well-behaved during the initial part of the de- 
cay. Furthermore, the statistical spread of the results 
(exemplified by the discrepancy between the median and 
the mean) is reduced as we approach the weakly-coupled 
bath regime. 

Thermal polarization of an electron spin bath is of- 
ten feasible with the temperatures and magnetic fields 
typically proposed in solid state spin quantum informa- 
tion processing. At the standard g = 2 value, the elec- 
tron Zeeman splitting corresponds to 1.3 K per Tesla. In 



FIG. 11: Comparison of Hahn echo results for different de- 
grees of thermal polarization: zero polarization (left), B /T — 
1 Tesla / Kelvin (center), and B/T = 2 Tesla / Kelvin 
(right). The zero polarization, infinite temperature results 
are shown in blue (color online), provided for compariza- 
tion along the polarization columns. The upper panels dis- 
play ^cCE (+' s )' ^CCE (triangles) and ^CCE (squares) mean 
value results. The lower panels display ^^CE me dian (di- 
amonds) as well as mean (squares) value results on a loga- 
rithmic scale. The £ scaling parameter is used as in Fig. [8] 
g 2 C E = C x 4 x 10 13 /cm 3 . 



fects of polarization are fairly straightforward. Since we 
are taking the limit of large applied magnetic field and 
using the secular approximation [see Eq. (|3|], the number 
of pairs that may contribute at the 2-cluster level is pro- 
portional to the probability of each spin being up times 
the probability of each spin being down: A^ pairs ex p^p± 
with p t/i = exp(±E z /2k B T)/2cosh(E z /2k B T) accord- 
ing to Boltzmann statistics where E z is the electron Zee- 

(2) 

man energy splitting. Since I/^ce * s a P r °duct over con- 

(2) 

tributing pairs, In (^cce) ^ PtPi- Note, however, that 
contributions from larger clusters play a more significant 
role earlier in the decay as we increase the polarization. 

Now we consider variations pertaining to resonances. 
In our canonical problem, all spins except for the central 
spin are taken to be resonant with each other. What 
happens when the central spin is resonant with the bath 
spins? In this case, we replace Eq. Q with 
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Fig. 11, we compare results for various polarized versions 
of our canonical problem. At the 2-cluster level, the ef- 



which allows flip-flops with the central spin. This induces 
both depolarization and dephasing together (i.e., where 
T2 may be T\ limited). These direct flip-flops produce 
non-trivial 1-cluster contributions in the CCE (depicted 
on the lower left of Fig. [l2| ). At short times, this dom- 
inates the spin echo decay. The spin echo that results 
in this scenario actually exhibit oscillations at the Zee- 
man precessional frequency (i.e., induced by the exter- 
nally applied B field). These oscillations are extremely 
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FIG. 12: (Top) Comparison of Hahn echo results for different 
resonance scenarios. We compare one bath species (left) and 
two bath species (right) with a 50/50 random mixture [see 
Eq. (14)]. Furthermore, results when the central spin is res- 
onant with bath spins [see Eq. (|13|)] are shown in blue. The 



upper panels display £qq E (filled circles) £qc E (+'s), and 
^CCE (squares) mean value results. Below those, the panels 
display ^^CE me dian (diamonds) as well as mean (squares) 
value results on a logarithmic scale. The £ scaling parameter 



(2) 



is used as in Fig. [8] g 2 C E = (x4x 10 13 /cm 3 . (Bottom) 
Depiction of a 1-cluster process when the bath spins are reso- 
nant with the central spin (left) and of an important type of 
4-cluster process when there are two bath species (right). 



fast, about 30 GHz at g = 2 and B = 1 Tesla. For 
numerical stability, we computed the CCE in this case 
by averaging over a few of these oscillation about each 
evaluated spin echo time; we define L$ in this manner. 

Relevant to many of our applications, another variant 
is to have multiple bath species that are only resonant 
within respective species. This applies to a bath of elec- 
trons bound to donor nuclei that have non-zero spin. For 
example, phosphorus nuclei have a 1/2 spin magnitude; 
this yields two bath species because up and down nuclei 
generate opposite hyperfine shifts for the electron spins. 
For this case, in the large applied magnetic field limit, 
we have 

ijeA 

-2j2KjS?S*. (14) 

Note that this is not quite the same as the combined inde- 
pendent effects of two baths at half concentration because 
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FIG. 13: (Top) Hahn echo results on logarithmic scales for 
a central electron spin at various distances from a sheet of 
random electron spins at a density of £ 2 ^ 3 x 10 n /cm 2 with a 
magnetic field parallel (left) and perpendicular (right) to the 
sheet. The distances are nm (black), C" 1/3 x20 nm 

(red), C~ 1/3 x 40 nm (green), and C~ 1/3 x 80 nm (blue), as 
labeled and generally from left to right (central spins more dis- 
tant from the surface have longer coherence times). The up- 
per panels display £qq E (+'s), ^cce (triangles), and £qq E 
(squares) mean value results. Below those, panels display 
^CCE median (diamonds) as well as mean (squares) value 
results. Bottom: Contour plots show relative strengths of 
dipolar couplings to the central spin from points on the sheet 
of bath spins (darker regions have larger relative coupling 
strengths) . 



there are inter-species Ising-like interactions. We com- 
pare various resonance scenarios (one versus two species 
and the central spin being resonant versus off-resonant 
with bath spins) in Fig. 12 along with schematic depic- 



tions of important process. 



B. Bath Geometry 

High concentrations of impurity spins may occur at 
the interface between materials. For example, dangling 
chemical bonds may host unpaired electrons. As a vari- 
ation of our canonical problem with this in mind, we 
consider two-dimensional geometries of bath spins. Our 
central spin may be at various distances (depths) from 
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this sheet of bath spins. We consider the limit of a large 
applied magnetic field, but the results will depend upon 
the angle of this applied field relative to the sheet of bath 



spins. We show Hahn echo results in Fig. 13 comparing 
various central spin depths for a magnetic field that is 
parallel or perpendicular to a sheet of bath spins. The 
spin echo curves in the parallel case are fairly smooth 
compared with the perpendicular case; this is due to dif- 
ferences in the spatial dependence of the dipolar coupling 
to the central spin for the different magnetic field angles 
(see the bottom of Fig. 13). 

Since our effective Hamiltonian [Eq. ([3|] scales in- 
versely with distance cubed (dipolar interactions), we use 
a scale factor £ for rescaling time, concentration, and 
depth together appropriately. £ = 1 is for a bath concen- 
tration of 10 n /cm 2 . These results are therefore applica- 
ble to various bath concentrations at the appropriately 
rescaled central spin depths. 



Finite spatial extent of the central spin electron 
wave function 



Our canonical problem idealizes the central spin and 
bath spins as being localized with zero extent (points) for 
the purposes of computing the dipolar interactions [see 
Eq. ([2])]. This is a reasonable approximation for a sparse 
system of donor-bound electrons such as Si:P. However, 
electrostatically defined quantum dots may have consid- 
erable lateral extent (e.g., roughly 50 nm in Ref . l46j) . To 
explore the impact of this finite extent of the central elec- 
tron's wave function, we use a simple Gaussian-shaped 
wave-function model in which the relative probability of 
electron occupation is given by 



P(x) oc exp 



cos 



(*!), (15, 



for \xs\ < 5/2 and P(x) = otherwise. We use 
ro = 50 nm as a defining lateral radius and 5 = 5 nm 
as a defining thickness. The coordinates are labeled with 
1, 2, and 3 so they will not be confused with the co- 
ordinate system of Eqs. ([I]) and ([3| which define z to 
point in the direction of the applied magnetic field. The 
xs direction is normal to an imagined surface forming a 
two-dimensional electron gas out of which the quantum 
dot is isolated. 

The finite extent of the electron's wave function im- 
pacts the dipolar interactions between the central spin 
and any bath spin. We may still use the effective Hamil- 
tonian of Eq. (J3j) , but must redefine the 6q,z = ^,0 such 
that 



-1 



(16) 



where is the position of the nth bath spin relative to 
the center of the quantum dot and we define 

(D atP ) (r) = J dxP(x)£> a>/3 (r - x). (17) 
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FIG. 14: (color online) Color map re pres entation of the 
non-zero dipolar tensor elements of Eq. (17) in the X2 = 
plane for interactions with a quantum dot electron geomet- 
rically defined by Eq. (15). By the x\ versus X2 symmetry, 
the X2 = plane is sufficient to express the dipolar infor- 
mation. Furthermore, due to the reflection symmetry about 
the x\ and X3 axes, one quadrant is sufficient for each of the 
non-zero tensor elements. We multiply the tensorial values 
by x 3 = (xl + xl) 3 ' 2 for convenience. The probability den- 
sity of the electron in the X2 = plane is represented in a 
grayscale image (black for high probability ranging to white 
for low probability) between the upper and lower quadrants. 



Here, D z z is an element of the dipolar tensor of Eq. Q 
where z refers to the direction of the applied magnetic 
field. For full generality, we compute the (D aj p) tensor 
in the quantum dot coordinate system and then rotate 
as appropriate to obtain (D ZfZ ). Due to the symmetry 
between x\ and X2 in our round quantum dot defined by 
Eq. (15), we may express all of the tensorial information 
in the X2 = plane. In this plane, it is clear from Eq. ([2| 
that (D h2 ) = (L> 2 ,i) = <£>2,3> = <£>3,2> = 0. The re- 
maining non-zero tensor elements, (^2,2), (^3,3)5 
and (Di^) = (^1,3) are displayed in color map form in 
Each tensor element is represented in a different 



Fig. 14 



quadrant of the x 2 — plane; due to reflection symmetry 
about the x\ axis and the £3 axis, one quadrant each is 
sufficient to convey the information for all quadrants. 

We generated the information that is represented in 
Fig. 14 from Eq. (17) using integration by Monte Carlo 



sampling at each point, r, on a two-dimensional grid of 
the X2 = plane. Using data generated in this fashion, 
and linear interpolation between grid points, we com- 
pute the Hahn spin echo of the quantum dot amongst 
various concentrations of a bath of point dipole electron 
spins in a 3-D bath (despite using a quantum dot typ- 
ically defined from 2DEGs). Results are shown in the 



upper panels of Fig. [15] for two different magnetic field 
directions: parallel to the surface (e.g., along x\ or #2), 
and perpendicular to the surface (i.e., along x%). To en- 
sure adequate grid spacing and precision for our dipolar 
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tensor data, we compared our Hahn echo results with re- 
sults using independent dipolar tensor data generated for 
larger grid spacing by a factor of two and found there to 
be no significant difference. 

In our canonical problem, we found that a rescaling of 
the concentration of bath spins is equivalent to an inverse 
rescaling of time. In Fig. [l5j we show the deviation from 
this behavior for our quantum dot qubit by plotting Hahn 
echo data versus a time scale that adjusts inversely with 
the concentration of bath spins. This deviation is not 
terribly large up to C E = 10 15 /cm 3 . At C E = 10 16 /cm 3 , 
the spatial extent of the central spin causes decoherence 
that has a significantly faster initial decay than a point 
dipole central spin. This can be understood by consider- 
ing the enhancement of dipolar interaction strengths for 
bath spins that are near some part of the laterally ex- 
tended quantum dot region but not particularly close to 
the center of this region. However, in some regimes, such 
as at later times, the lateral extent actually causes en- 
hanced coherence. This counterintuitive enhancement is 
an effect of the anisotropy of the dipolar interactions. It 
is erratic and differs for the two cases of the differing In- 
field directions. The lower panel of Fig. 15 shows results 
using an artificial 1/R? potential for the interactions with 
the central spin, removing the anisotropy, and we find a 
consistent decrease in the relative coherence, quantum 
dot central spin versus point dipole central spin, as the 
bath concentration is increased. Thus, the intuitive un- 
derstanding is only thwarted, in some regimes and only 
slightly, by the anisotropy of the dipolar interactions. 



B-field Parallel to Surface B-field Perpendicular to Surface 



D. 



Bath geometry and the wave function of the 
central spin 



Combining bath geometry considerations of Sec. IV B 
with finite extent of the central electron's wave function 
discussed in Sec. IV C[ we see a variety of different trends. 
In Fig. [l6j we compare the decoherence effects from a 
two-dimensional sheet of bath spins at random positions 
for a point-like central spin and a Gaussian-shaped quan- 
tum dot central spin with the bath sheet at a distance of 
3 nm from the center of the central spin. For the quan- 
tum dot central spin case, this corresponds to 5 Aabove 
the edge of the wave function, xs = 5/2 + 5 A[Eq. (fl5|]; 
we are considering this roughly as a limiting case for a 
sheet of electron spins that is very close to the quantum 
dot but spatially independent. Unlike the trend observed 
when changing the qubit's wave function in the canon- 
ical three-dimensional sparse bath (Fig. 15), the coher- 
ence time for the quantum dot central spin is actually 
increased relative to the point-like central spin. In the 
former case, more decoherence resulted from enhanced 
dipolar interactions near the laterally extended quantum 
dot region. Here, where we consider a 2-D bath at densi- 
ties of 10 11 /cm 2 or 10 12 /cm 2 , the extent of the quantum 
dot actually reduced decoherence because it tends to re- 
duce its sensitivity to flip-flopping pairs of nearby bath 
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FIG. 15: (Top) Comparison of Hahn echo results affected by 
the spatial extent of the central spin electron to various de- 
grees. Given a central spin "quantum dot" geometry defined 
by Eq. (15), we show results for various electron spin bath 



concentrations, Ce, and compare with the canonical results in 
which the central spin has zero extent. The applied magnetic 
field is parallel (left) or perpendicular (r ight ) to the surface 
whose normal is the X3 direction of Eq. (15). We use g = 2 
(g 2 CE — C x 4 x 10 13 /cm 3 ). Using the scaling parameter £, 
time is rescaled inversely with the bath concentration so that 
the "zero extent" curves are universal for all concentrations. 
The upper panels display £qq E (+'s), and £qq E (squares) 
mean value results. Below those, panels display ^^CE median 
(dashed line with diamonds) as well as mean (solid line with 
squares) value results. (Bottom) Corresponding results when 
we use an artificial isotropic 1/R 3 interaction with the central 
spin. £qq E mean value (solid +'s) and median value (dashed 
x's) results are displayed. Showing only 2-cluster results is 
sufficient to make our point: without the anisotropy of the 
interactions there is a consistent trend, monotonic in £. 
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FIG. 16: (Top) Hahn echo result of a Gaussian-shaped quan- 
tum dot defined by Eq. |l5| with decoherence induced from a 
two-dimensional bath, in the xi,X2 plane, 5 A from its edge 
(xs = 3 nm). This was chosen as somewhat of a limiting 
case in terms of the closeness of the bath to the central spin 
quantum dot. Dashed lines show results for a point-like cen- 
tral spin for comparison. We show mean value L^ CE (+'s), 
^CCE (triangles), and ^^CE ( sc t uares ) results, as well as me- 
dian value ^cce (diamonds) results. Different colors (color 
online) are used to help make the curves more easily distin- 
guishable. Bath densities are 10 n /cm 2 (left) and 10 12 /cm 2 
(right), and the applied magnetic field is perpendicular to the 
bath and dot (upper) or parallel to the bath and dot (lower). 
(Bottom) Depiction of a laterally extended quantum dot in 
the presence of sheet of bath spins. 
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FIG. 17: Hahn ech o re sult of a Gaussian-shaped quantum 
dot defined by Eq. (15) with decoherence induced from a 



two-dimensional bath, in the x±,X2 plane, at various depth 
distances along X3. Dotted lines show results for a zero ex- 
tent (point-like) central spin for comparison. Bath densities 
are 10 n /cm 2 (left) and 10 12 /cm 2 (right), and the applied 
magnetic field is perpendicular to the bath and dot (top) 
or parallel to the bath and dot (bottom). As labeled, the 



depths (color online) are defined in correspondence to Fig. [13 
1/3 x 10 nm (black), C 



c 



-1/3 



x 20 nm (red), C" 1/3 x 40 nm 



(green), ( 1/3 x 80 nm (blue), ( 1/3 x 160 nm (orange), and 



" 1/3 x 320 nm (purple) where £ = 1 for 10 n /cm 2 and 
C 10 3/2 « 32 for 10 12 /cm 2 (C~ 1/3 10~ 1/2 « 0.32). Ad- 
ditionally, the 3 nm depth results from Fig. [16] are displayed 
in thick brown. These are all results of mean values 

(+'s) as a reasonable short time approximation to understand 
trends. 



spins. To put this another way, the important factor is 
the difference in the effective magnetic field experienced 
by the central spin as bath spins flip-flop, not the absolute 
magnitude of their interactions. Thus, while the extent 
of the quantum dot wave function tends to increase the 
strength of dipolar interactions to the bath spins, in cer- 
tain geometries, this difference in interactions amongst 
nearby bath spins may be reduced. The depiction at 
the bottom of Fig. [16] helps to illustrate this effect. It 



should also be noted that, in going from the point-like 
limit to the extended quantum dot wave function, we en- 
ter a regime in which the mean and median of the Hahn 
echo decay are nearly identical; indeed, for a sufficiently 
dense bath, we expect the decay to result from a large 
number of small cluster contributions to govern the de- 
coherence and, as an effect of the central limit theorem, 
most random instances of the bath should and do yield 
similar results. 

In Fig. [lTj we show decoherence induced by 2-D sheets 
of randomly located bath spins at various depths from a 
quantum dot central spin. We show only the 2-cluster 



results as a rough approximation to understand general 
trends as we approach the large distance limit in which 
the zero-extent approximation of the central spin is valid. 
Some peculiarities emerge. First, in all of the cases 
we consider [10 n /cm 2 (left) and 10 12 /cm 2 and different 
magnetic field orientations], there is actually an initial 
decrease in coherence time as we move beyond a 3 nm 
depth from the bath. This is somewhat counterintuitive, 
but the effects of the change in dipolar interactions with 
changing depth is not so straightforward due to the angu- 
lar dependence of the dipolar interactions. Beyond about 
10 nm, however, it follows the more intuitive trend that 
increasing the depth causes coherence times to increase. 
There is an exception for the case with a 10 12 /cm 2 bath 
density and parallel B field in going from C" 1/3 x 160 nm 
to C _1//3 x 320 nm depth (51 nm to 101 nm); this is a 
relatively minor effect that is apparently related to the 
anisotropy of dipolar interactions. Generally we find that 
we approach the large distance limit at a depth of about 
C" 1/3 x 320 nm for the two densities we study. In this 
limit, we can neglect the spatial extent of the central spin 
wave function and use the results from Sec. IIVBI 
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E. Pulse Sequences 



As we discussed in Sec. |IIB[ we chose an ideal Hahn 
spin echo as the context of our canonical problem as a 
standard approach of eliminating the effects of inhomo- 
geneous broadening. In this section, we consider other 
pulse sequence scenarios. We first examine the effects 
of inhomogeneous broadening itself which can incur er- 
rors when refocusing pulses are not used. We display the 
probability distribution of magnetic field shifts of an en- 
semble, responsible for inhomogeneous broadening, due 
to background electron spins in the infinite temperature, 
point dipole limit on the upper left side of Fig. [18} On 



the upper right side of Fig.[18j we compare coherence ver- 
sus time calculations in four different contexts: ensemble 
free induction decay (FID), FID with characterization, 
Hahn spin ec ho, an d a two-pulse Carr-Purcell-Meiboom- 
Gill (CPMGjSEH sequence. 

The ensemble FID corresponds to the previously de- 
scribed case of ensemble averaging of the signal which 
has no dependence upon bath dynamics whatsoever. FID 
with characterization, or narrowed state FID (NFID), 
takes the inhomogeneity out of the problem by account- 
ing for a static offset in the magnetic field character- 
ized individually for each qubit; coherence of the central 
spin is lost, however, due to dynamics of the bathP^EU 
This corresponds to an experiment in which the central 
spin splitting is pre-measured, or prepared with a well 
fixed value. The Hahn spin echo is simply the case of 
the canonical problem. In the context of ideal applied 
pulses in a dephasing-only limit, the CPMG is equiva- 
lent to the Carr-Purcell sequence and is represented as 
(r 7T r) n with t = 2nr. The two-pulse CPMG se- 
quence, t— )>7r— )>2t— )>7r— ^r, is equivalent to the two- 
pulse Uhrig dynamical decouplings 5 -^ (UDD) sequence 
as well as a first-level concatenation of the spin echo se- 
quenced^ Our CCE methods, as discussed in Sec. Ill 
work well for FID with characterization and for the two- 
pulse CPMG down to below 50% of the decay. For these 
cases, we do not show the tail of the decay for Lqq E mean 
results because we fail to obtain convergent results. How- 
ever, the initial part of the decay is often most relevant 
for quantum computation in any case. 

The coherence time improves successively in the FID- 
NFID-SE-CPMG 2 sequence of experiments. In Refs. [52] 
and l42| CDD sequences were shown to work well for spec- 
tral diffusion in cases where the intra-bath coupling is 
weak compared with the qubit-bath coupling and UDD 
sequences work even better in cases where the relevant 
time scale is short compared with all of the couplings. 
Neither of these perturbations are particularly relevant 
in the problem considered here. If we were in the regime 
of weak intra-bath coupling, larger clusters (than pairs) 
might dominate the two-pulse CPMG results as observed 
in Refs. l52l and [42l This is not observed here. If we were 
in the short time regime, we would see an exp (— t 4 ) de- 
cay for the Hahn spin echo and exp (— t 6 ) decay for the 
two-pulse CPMG. We do not; in the logarithmic-scale 



plots of Fig. [THJ we do not observe a significant change 
of slope between the Hahn spin echo decay and the two- 
pulse CPMG decay. We do see some improvement with 
the two-pulse CPMG sequence in any case; although, as 
a function of the time between pulses r, we roughly break 
even relative to the Hahn spin echo (bottom of Fig. 18). 
This is consistent with the calculations based on classi- 
cal Ornstein-Uhlenbeck noise: the 1/uj 2 tail of the spec- 
tral density of this noise leads to exp(— t 3 ) decay for any 
dynamical decoupling pulse sequence, and it leads to a 
sublinear T2 ~ n 2 / 3 s caling of the coherence time with 
the number of pulsesPSEME] 



V. APPLICATIONS 

Now that we have presented a detailed study of our 
canonical problem in Sees. |H| and |III| an d looked at sev- 
eral variants of the problem in Sec. |IV[ we now discuss 
specific applications. We consider two crystalline mate- 
rial substrates, silicon and carbon (in diamond form), in 
which nuclear spins may be nearly eliminated through 
enrichment. In both of these cases, however, impurity 
electron spins can become the predominant source of de- 
coherence. Such spin baths relate to our canonical prob- 
lem and its variants. 



A. Donor in silicon 

Donors in silicon make a promising candidates to host 
quantum bits in the form of electron spins. In bulk, ex- 
periments indicate that donor-bound electrons can main- 
tain spin coherence on the timescale of a second by us- 
ing enriched silicon with few nuclear spins and single- 
spin read-outP^ as well as coherent controP^"^ has been 
demonstrated. In this section, we consider common spin 
baths that may affect these qubits: 29 Si, background 
phosphorus dono rs, and electron spins at an interface. 

Previous workP^^ has demonstrated remarkable 
agreement of cluster expansion calculations for spin echo 
decay with corresponding elect ron spin resonant mea- 
surements in natural silicon^ 2 -^ and accurately predicted 
the effects of varyi ng 29 Si concentrations with experimen- 
tal measurements, ^ 14 * 15 * . In Ref. 23, we examined the 
effects of background phosphorus donors when the 29 Si 
has been reduced to very low concentrations through en- 
richment and demonstrated agreement with experiments 
now published in Ref. |25j 

Ignoring the background phosphorus donors for a mo- 
ment, we note that a spin bath of 29 Si in the low concen- 
tration limit, in which the spins act as point dipoles, cor- 



responds to the <^>o < go variant of Fig. 10 In this case, 
go ~ 3000 x^>o; the electron spin has a much larger mag- 
netic moment than the 29 Si. The background phosphorus 
donor spin bath alone, in the low concentration limit, 
corresponds to the two bath species variant of Fig. 12 
because of the spin-1/2 phosphorus nucleus; to a good 
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FIG. 18: Top left: probability distribution of magnetic field 
shifts due to background electron spins in the infinite tem- 
perature, point dipole limit. The distribution is over differ- 
ent random spatial realizations and spin states of the back- 
ground spins. Top right: coherence versus time for ensemble 
FID (black) , FID with characterization (red) , Hahn spin echo 
(green), and CPMG (blue), with a left to right trend, re- 
spectively. We show mean values of £qq E (solid +'s) and 

^CCE ( SOU( l squares). The bottom of the right panels dis- 
plays results on a logarithmic scale and includes median val- 
ues of £qq E (dashed diamonds). Bottom: coherence versus 
r, the time between tt pulse for the same pulse sequences. 
r — t/2 for the Hahn spin echo and r — t/4 for the two- 
pulse CPMG. The £ scaling parameter is used as in Fig. [8] 
g 2 C E - C x 4 x 10 13 /cm 3 . 



approximation, we can assume that donors with differ- 
ing nuclear polarizations do not flip-flop with each other. 
In Fig. [l9j we show Hahn spin echo T2 times as a func- 
tion of ^3i (nuclear spin) concentration, Cat, in combi- 
nation with various concentrations of phosphorus donors 
(electron spins), Ce- The 29 Si bath induces decoher- 



29 Si concentration, C N (cm 3 ) 




FIG. 19: Hahn spin echo T2 times, when the Hahn echo 
reaches a value of exp (—1), for phosphorus donors as a func- 
tion of the fraction of 29 Si donors, or corresponding concentra- 
tion Cat, as well as the concentration of background donors 
Ce- At high Cat, contact hyperfine interactions dominate 
and T2 is dependent upon the magnetic field direction rela- 
tive to the lattice orientation. At low Cn, T2 is dependent 
upon Ce, and eventually dominated only by dipolar inter- 
actions (which includes dipolar-approximated electronuclear 
interactions). Experimental results are shown as square sym- 
bols, from Ref. 1151 and a star symbol at 50 ppm 29 Si, from 
Ref. 1251 This figure is slightly revised from that of Ref. 1231 
updating the 50 ppm 29 Si experimental value to the published 
T 2 450 ms. 



ence through its flip-flopping dynamics, but it also sup- 
presses donor-induced decoherence via Overhauser shifts 
that cause the donors to be off-resonant with each other. 
This effect is demonstrated in Fig. [19] by the initial in- 
crease of T2 as 29 Si is increased for the Ce > cases. 
However, this effect is not very prominent in the short 
time regime that is of most interest for meeting quantum 
error correction thresholds! 2 ^ 

In Fig.[20j we examine the Hahn spin echo correspond- 
ing with the T2 ~ 450 ms experiment from Ref. |25] This 
is a bulk experiment in which all the donors are sub- 
ject to the spin echo refocusing pulses. As a result, the 
measured decoherence signal is dominated by instanta- 
neous diffusion, the inhomogeneous broadening effects 
of the donors that are mutually flipped by the refocus- 
ing pulse. In other words, the spin echo is only able to 
cancel the effects of inhomogeneous broadening from the 
parts of the bath that remain unchanged by the refo- 
cusing pulse. However, by performing a series of mea- 
surements in which the angle of the refocusing pulse is 
varied (i.e., smaller than ir pulses), they extrapolate the 
effective single- spin decoherence. Our calculations ap- 
ply to this extrapolated result. In order to confirm that 
the dominant decoherence is due to flip-flopping donors, 
Ref. [25] presents measurements of donor decoherence in 
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FIG. 20: Hahn spin echo versus time corresponding with the 
50 ppm 29 Si and C E = 1.2 x 10 14 /cm 3 experiment of Ref. [25] 
Time is shown on a logarithmic scale in the main plot and a 
linear scale in the inset. The theory exhibits a T2 ~ 600 ms, 
reasonably close to the 450 ms T2 fit of the experiment! 2 -^. 
The short-time (first few percent of the decay) behavior is 
dominated by 1-cluster contributions shown with filled circle 
symbols; these are direct flip-flops between the central spin 
and other donors. Indirect flip-flop processes then dominate 
the decay, 2-clusters (+'s) then 4-clusters (squares). In the 
T 2 regime near 600 ms, the decoherence is a combined effect 
from 29 Si flip-flops as well as donor flip-flops. 



the presence of a significant magnetic field gradient, pro- 
longing T2 to an extrapolated value of about 10 s. Such a 
gradient suppresses donor flip-flopping by shifting them 
off resonance from each other. 

Two types of decoherence induced by donor flip- 
flopping are distinguished in Ref. [25] as indirect flip-flops 
(flip-flops among bath spins) and direct flip-flops (flip- 
flops with the central spin). In Ref. [23] we did not con- 
sider the effects of direct flip-flops between the central 
spin and background donor spins. This effect, which we 
examine in Fig. [l2] (where the central spin is resonant 
with bath spins), results in 1-cluster contributions that 
dominate in the short time regime. We include both di- 



rect and indirect flip-flops in Fig. [20] and note that the 
direct flip-flops have a negligible effect on the T2 time but 
significant effect on the short time behavior. 

It is important to keep in mind that there will be large 
statistical variations for different spatial realizations of 
spin baths when we are in the low concentration regimeP^ 
This has been indicated by differences in the median and 
mean values throughout this paper. The statistics are 
examined in greater detail in Fig. 21 for 29 Si spin baths 
and for background donor spin baths separately. 

In addition to 29 Si and background donors, interfaces 
that play important roles in semiconductor technology 
may introduces additional baths of spins that may in- 
duce decoherence. Experiments 60 have demonstrated 
that donors closest to an interface have shorter coher- 



FIG. 21: Statistical dependence of the Hahn spin echo upon 
various random instances of spin baths, 29 Si (left) or donors 
(right), in the low concentration limit. Top: Hahn echo er- 
ror (one minus the echo) for various percentiles, computed 
at each time point independently, and the ensemble averaged 
spin echo. Bottom: fits of the error to 1 — exp [— (t/T q ) nq ] « 
(t/T q ) ng in the 10 -4 error regime (motivated by common 
fault-tolerance thresholds) for various percentiles. The £ scal- 
ing parameter applies simultaneously to the time scale and 
concentration scale, having an inverse relationship in the low 
concentration limit. 



ence times which have thus far proven to be much shorter 
than the coherence times observed in the bulk. We re- 
fer to Fig. [13] with spin echo calculations for decoherence 
induced by a 2-D bath of electron spins. This is ap- 
plicable where we can approximate spin interactions as 
point dipole interactions and where bath spins are reso- 
nant with one another. Dangling bond spins at an inter- 
face may have ^-factor variations that cause them to be 
off-resonant with one another and suppress flip-flopping 
noise. The experiments of Ref. 60 exhibit faster spin echo 
decay than we would expect from calculations along the 
lines of Fig. [l3j even assuming no ^-factor variation. A 
different theory of dangling bond spin -induced decoher- 
ence proposed in Ref. 61, involving phonons and spin- 
orbit interactions, matches with Ref. [60] but it requires 
dangling bond concentrations as high as 10 13 /cm 2 , which 
might be unrealistic. There is no consensus at this time 
as to the true cause of the decoherence in Ref. 60. One 
suspicion is that the noise is induced by exchange-coupled 
parasitic dots at the interface such as observed in Ref. [62] 



B. Nuclear spin qubit 

It has been proposed^ to use donor nuclei for quan- 
tum memory storage in between quantum information 
processing. Nuclear spins have much weaker magnetic 
moments than electrons and are thus less susceptible to 
magnetic field noise, leading to longer coherence times. 
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FIG. 22: ^cCE Hahn spin echo results for a phosphorus 
nuclear spin in a bath of 29 Si. Mean values are encircled +'s 
and median values are encircled x's. The left panels are for 
natural Si with 4.67% 29 Si and includes lattice effects. The 
right panels are for varied low concentrations of 29 Si in the 
continuum limit. The bottom panels present the same data 
as respective top panels but in a logarithmic scale. 



It is interesting to note that the scenario of a central nu- 
clear spin in a nuclear spin bath, where the interactions 
strengths with the central spin and among the bath are 
comparable, is essentially the same as a central electron 
spin in an electron spin bath. Thus, our methods are 
applicable to nuclear spin decoherence that is caused by 
other nuclei. In Fig. [22] we present calculated results 
for the Hahn echo decay of a nuclear spin of a phos- 
phorus donor in silicon with 29 Si as bath spins. The 
left panels are for natural Si and include lattice effects 
while the right panels apply to varied 29 Si concentra- 
tions in the continuum limit. The range of values for 
different spatial configurations of the 29 Si are consistent 
with recently reported single P nucleus T 2 measurement 
of roughly 60 ms (see Ref.[58]) and 30 ms (see Refs.l64land 
[65]) . The right panels are analogous to the 2-cluster part 
of Fig. [8j however, it differs from the canonical problem 
in that the central spin gyromagnetic ratio is roughly a 
factor of two larger than that of the bath spins. The vari- 
ant that is explored in Fig. [l0]of Sec. |IV A| is applicable 
here. For that variant, we found that the cluster expan- 
sion has better convergence when the central spin has a 
larger gyromagnetic ratio (e.g., g-factor) than the bath 
spins. For that reason, the decay is well approximated 
with 2-clusters. 



C. Quantum dot in silicon 

Quantum dot electron spins are promising for qubit re- 
alizations. In particular, using a singlet /triplet encoding 
for two electrons in double quantum dots, demonstrated 
in GaAs^^^, fast single-qubit operations with elec- 



ppm 



FIG. 23: Computed Hahn spin echo T2 times, when the Hahn 
echo reaches a values of exp ( — 1), as well as T 2 * times, when 
the ensemble FID decay reaches exp (—1), for P donors as well 
as qua ntum dots (q-dots) defined with the probability density 
of Eq. ( 15 ) in Si. We consider decoherence induced by various 
concentrations of 29 Si. For our quantum dot model, we chose 
the lattice orientation such that [100] corresponds with the 
X3 direction (normal to the confinement well) as depicted. 



trical controls and readoutP^^ are possible. In silicon, 
quantum dot spin qubits are starting to be realizedp3E3 
and long coherence times are possible, particularly with 
isotopic enrichment of Si^ (and also Ge in Si/SiGe quan- 
tum dots^). 

Even with very high isotopic enrichment, as with the 
silicon donor qubits discussed in Sec. V A , coherence 
times will be limited by a background of impurity phos- 
phorus donors. In Sec. |IVC| we found that the effects of 
the spatial extent of a quantum dot electron in such a spin 
bath are negligible up to concentrations of 10 14 /cm 3 . So 
at very achievable levels of silicon purity, the wave func- 
tion extent is insigificant and the background phosphorus 
induced dec ohere nce problem is essentially the same as 
that of Sec. VA and the two-species variant (due to the 



two P donor nuclear polarizations) of the canonical prob- 
lem as presented in Fig. [12] 

In addition to background phosphorus donors, 2-D 
electron spin baths may be present. For example, 
Si/SiGe quantum dot structure may employ modulation 
doping layers with some fraction of un-ionized donors. 
For these effects, we refer to Sec. |IVD| where we study 
electron spin decoherence for a quantum dot extended 
wave function with a 2-D electron spin bath at various 
distances. This study has important implications for tol- 
erated densities and distances of such a bath in order to 
achieve desired decoherence times. 

It is also important to understand how the effects of 
various 29 Si concentrations may differ for a quantum dot 
compared with the donor qubit discussed in Sec. VA In 
Fig. 



23 , we present decoherence time versus 29 Si concen 
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tration in contrast with results for the donor qubits and 
observe the effects of the difference in the qubit wave- 
function shape. At very low concentrations, the wave- 
function shape is irrelevant and a point-like model works 
for either type of qubit. At moderate concentrations, the 
laterally extended quantum dot qubit has increased cou- 
pling to a number of bath spins so that decoherence times 
are faster than those of the donor qubit. This is the same 
effect observed in Fig.[l4]of Sec. |IV C] for an extended cen- 
tral electron spin in an electron spin bath at 10 15 /cm 3 . 
At higher concentrations, the situation is reversed. The 
effect here is analogous to what was observed in Fig. 16 



2-cluster only 



up to 4-clusters 



of Sec. IV D for the effects of a laterally extended central 
spin near a 2-D bath; the common cause is a reduced 
sensitivity to flip-flopping of neighboring bath spins. 

At higher concentration, the decoherence is affected 
by the magnetic field angle relative to the lattice orien- 
tation, having the shortest coherence times when the B 
field is aligned with the nearest neighbor direction. For 
the quantum dot wave function, this effect is less pro- 
nouned compared with P donors but is still appreciable. 
Around 2 to 4 ppm of 29 Si, the quantum dot case with B 
parallel to [111] is affected by the anisotropy of the dipo- 
lar interactions in an unusual way that causes T2 to be 
long compared with other 5-field directions. With just 
these few exceptions, however, T2 is not strongly depen- 
dent upon the 5-field direction. 



We also plot T| in Fig. 23 for P donors and quan- 
tum dots. Above roughly 30 ppm 29 Si, the quantum dot 
exhibits longer T 2 * times than the P donors as expected 
from the central limit theorem that predicts longer T| for 
wave functions with greater spatial extent: T| oc y/~N. 
The central limit theorem, however, does not apply well 
to the donor case, particularly at low densities. In fact, 
the magnetic field shift probability distribution for point 
dipoles in the top left of Fig. [18] fits a Lorenzian distribu- 
tion much better than a Gaussian distribution. In the low 
concentration regime, the roles are reversed; T| is shorter 
for the quantum dot that has a few strong 29 Si contact 
hyperfine interactions than for the donor with negligible 
29 Si contact hyperfine interactions. At low enough con- 
centrations, the extent of the qubit wave function should 
be negligible in either case and the T 2 * times should be 
the same, but that regime is well below the range we 
present. It is interesting to note that we do approach 
that regime for T2 but not T 2 *. To understand this, con- 
sider that all bath spin give direct contributions to T 2 * 
while, in our large 5-field limit, bath spins contribute to 
T2 only indirectly via near-resonant flip-flopping pairs. 

The effect of the magnetic field angle is negligible for 
the T 2 * data we present. This is expected where the 
isotropic contact hyperfine interaction dominates. It is 
also expected in the low concentration limit where the 
extent of the wave function has a negligible upon T 2 *; 
in the presented range, this limit is approached for the 
donor but not the quantum dot. 
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FIG. 24: Analogous to Fig. [8] but with the spin echo results 

for an NV center in a bath of PI centers at a concentration 

of C E = C x 10 19 /cm 3 . Each "spaghetti" strand is the result 

for a different random spatial configuration averaged over a 

(2) 

large number of spin state templates. The mean of ^CE an< ^ 
^CCE are re P resen ted as encircled +'s and squares, respec- 
tively. The median of ^^CE and £qq E are represented as 
encircled x's, and diamonds, respectively. The bottom figures 
show the same spin echo results on a logarithmic scale for the 
decay. 



D. NV Center in Diamond 

Nitrogen vacancy (NV) defect centers in diamond form 
remarkable qubits which may b e coherently controlled at 
room temperature! 8 * 26 ^ 29 * 31 * 80 ! ^ Electron spins of nitro- 
gen atoms known as PI centers in typical diamond are 
a major source of decoherence. At concentrations of be- 
low 200 ppm (3.5 x 10 19 /cm 3 ), spin echo decoherence 
times of about T2 ~ 3/is have been observed at room 
temperature! 2 ^ High-purity diamond with low concentra- 
tions of nitrogen decohere due to 13 C nuclear spins with 
spin echo times of T 2 ~ 13 /isP^In isotopically enriched 
high-purity diamond with O.c 
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C and paramagnetic 
defects below 10 13 /cm 3 , spin echo coherence as long as 
T2 = 1.8 ms has been measured. 28 

The NV decoherence problem in the presence of PI 
centers in a large magnetic field is similar to our canoni- 
cal problem. There are a few important differences. De- 
tails of the interactions of the NV and PI center system 
appear in the supplemental online material of Ref. [27j 
The important features for our consideration are the fol- 
lowing. The NV center is treated as a localized spin 
with So = 1 and g = 2 (the free electron g-factor). The 
two-level qubit system in Ref. |27l and in our calcula- 
tion, is the m s = {0, —1} subspace. The PI centers are 
spin 1/2 electrons (with g = 2) bound to spin 1 nitro- 
gen nuclear spins. Each center has a derealization axis 
indicating the neighboring carbon atom is sharing the 
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electron with the nitrogen atom. This derealization axis 
changes over a time that is much longer than the charac- 
teristic time of a single experimental run in Ref. |27l In 
the large magnetic field regime, e.g., B — 740 G, as in 
Ref. [27j the derealization axis simply determines a hy- 
perfine energy coupling of the PI centers: A\ = 114 MHz 
for the [111] axis and A 1 = 86 MHz for [111], [ill], or 
[111] where AiS%Ij: is the hyperfine energy shift with 
Sk as a PI electron spin operator and 1^ as its nitrogen 
spin operator. We therefore have five different species 
in our spin bath corresponding to five different hyperfine 
shifts with various population percentages: 1/12 fraction 
each for Ai = 114 MHz and If = ±1, 1/4 fraction each 
for A x = 86 MHz and l z k = ±1, and 1/3 fraction for 
1^ = 0. For our high-field limit calculation, we neglect 
flip-flopping between spins of different species; these are 
off-resonant with each other. We can also safely neglect 
flip-flopping between the NV center and any bath spin. 
Our results in Fig. [24] show T2 ~ 3 /is for a PI cen- 
ter concentration of 10^ 9 /cm 3 , consistent with measure- 
ments reported in Refs. 1271 and 1291 where 10 19 — 10 20 /cm 3 
concentrations are estimated. 



VI. CONCLUSION 

We have presented a detailed account of a cluster-based 
theory of spin echo decoherence of a central spin (qubit) 
interacting dipolarly with a bath of spins of the same 
kind (i.e. we focused on the case of the qubit-bath cou- 
pling being the same as the intrabath coupling). While 
the previously developed cluster theories were proven to 
be very successful in the case of the intrabath coupling 
being much weaker than the qubit-bath coupling, the all- 
dipolar problem with symmetric couplings requires the 
cluster-based approach to be nontrivially modified. We 
have shown that the decay of the spin echo signal can 
be calculated reliably by solving for evolution of finite 
groups (clusters) of bath spins coupled to the central 
spin, provided that the offsets of the splittings of these 
spins cause by dipolar interactions with the rest of the 
bath are properly taken into account. In a sparse bath 
the disorder in energy splittings of bath spins leads to 
localization of flip-flop dynamics, i.e. at the timescale at 
which the qubit's coherence decays it is enough to con- 
sider the dynamics of still rather small clusters (up to 
4-6 spins). This result does not follow from any kind of 
simple perturbative argument, and while it could have 
been suspected, the existence of such an effect had to be 
checked by careful numerical simulations. 

Our theory allows for quantitative evaluation of de- 
coherence in an all-dipolar system of spins. It presents 
a microscopic (i.e. derived from the Hamiltonian of the 
system) solution to the original spectral diffusion prob- 
lem, which has been approached by phenomenological or 
semi-phenomenological stochastic theories for more than 
fifty years. We have presented a broad selection of re- 
alistic applications of this theory, including calculations 



of spin echo decay for (1) electrons bound to phospho- 
rus donors in isotopically purified silicon (reported pre- 
viously in Ref. 23 and recently confirmed experimentally 
in Ref. [25]) : (2) nuclear spin qubits in silicon; (3) quan- 
tum dots in isotopically purified silicon; and (4) nitrogen- 
vacancy (NV) centers in diamond, in the case in which 
the nitrogen spins are the dominant source of decoher- 
ence. Although we have mainly emphasized Hahn spin 
echo decoherence, our technique is very general and can 
be applied to any quantum control context in principle. 
Together with previous works on cluster theories of de- 
coherence due to dipolar interactions among the bath 
spins^nuJIHl ( m which the qubit-bath hyperfine coupling 
was much stronger than the intra-bath interaction), the 
theory from this paper completes a body of work de- 
voted to realistic calculations of decoherence in systems 
in which the dipolar interactions within the bath play a 
dominant role. 

The temporal nature of these decoherence problems is 
an important aspect for allowing cluster expansions to 
succeed. The perturbative arguments always have a fac- 
tor of time accompanying interaction energies. As time 
increases, the effective perturbation parameter increases, 
reducing the performance of the perturbative (cluster) 
expansion. At long times, we observe that the expansion 
fails entirely. We consider a cluster expansion successful, 
however, when it is well-behaved and convergent on the 
timescale of the decay. Furthermore, only the initial part 
of the decay is of interest for typical quantum computing 
applications. Although we do not prove formal conver- 
gence in the cluster expansion or our selection heuristics, 
we demonstrate good convergence in practice going out 
to the 6-cluster order of the expansion provided that sim- 
ulation times are sufficiently short. 

While a standard desktop computer is capable of pro- 
ducing many of the results that we present here in a rea- 
sonable amount of time (hours for a 2-cluster spin echo 
decay of a typical scenario with good accuracy), we made 
significant use of Sandia's high performance resources to 
acquire accurate results for the wide range of scenarios 
in our study as well as a lot of experimentation with 
our methods and heuristics. Runtimes increase signifi- 
cantly with increasing cluster size. The cluster selection 
heuristics we use are critical for making larger cluster cal- 
culations feasible. Fortunately, the calculations are easy 
to parallelize with each processor treating a particular 
instance of a random bath instantiation. 

Possible future directions will be to probe spectral den- 
sities of the bath directly using cluster techniques and 
to answer questions regarding conditions under which 
the bath may be treated independently from qubit con- 
trol (i.e., a classical bath). Spectral density descriptions 
are very convenient and are commonly probed in exper- 
iment. Whether or not a bath acts classically has im- 
portant implications for quantum control. The connec- 
tion to Ornstein-Uhlenbeck noise in our all-dipolar model 
that others have noted and we confirm does imply that a 
classical noise model is applicable in some manner (i.e., 
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for a particular spatial configuration averaged over spin 
states). 
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Appendix A: Cluster Sampling Heuristics 

The relative importance of clusters may be judged, 
heuristically, by the following two factors: the strength of 
coupling between the central and bath spins of the clus- 
ter, and the strength of coupling among the bath spins of 
the cluster. We address the first factor with a radial cut- 
off. Since interactions decrease with increasing distance, 
we'll ignore clusters that are fully outside of some cutoff 
radius, Rc relative to the central spin. For the remaining 
clusters, those with at least one bath spin within Rc from 
the central spin, we apply a heuristic to address the sec- 
ond consideration (the coupling strength among spins). 
Motivated by Lemma 1 of Sec. |IIIB[ which is not com- 
pletely valid when using the interlaced spin state averag- 



ing of Sec. HID but should still have some approximate 
validity, we want clusters that may be fully connected 
with sufficiently strong interactions. We will therefore 
assign a heuristic strength of a cluster to be the small- 
est interaction neccessary to complete the connectivity 
of the cluster (as in Fig. [2j. Algorithmically, we may 
compute this by successively picking off the strongest in- 
teraction among bath cluster spins until the cluster has 
full connectivity from the picked interactions; the last 
interaction picked (the smallest necessary interaction) is 
the heuristic strength of the cluster. 

In addition to Rc, we also employ a cutoff for the 
maximum number, TV^, of clusters to be considered for 
each cluster size, k. For each CCE computation (which is 
typically averaged over bath spin locations and initial po- 
larizations) and for each cluster size k (up to some max- 
imum), we select the clusters of the highest heuristic 
strength for evaluation. In the donor spectral diffusion 
application, we use one more cutoff in which we only 
consider, for the sake of cluster selection, intereractions 
among donors with comparable Overhauser shifts; this is 
a resonance energy cutoff, Ec- 



FIG. 25: Illustration of the constructive proof for Lemma 
2, showing a path, with numbered edges and vertices, that 
traverses the entire graph Q. A connected sub-graph 7i is 
formed from the edges and vertices of the path with the ex- 
clusion of the final edge and vertex. The edges and vertex 
that are excluded from 7i are in gray. 



The key to an efficient algorithm for finding each set of 
Nk clusters, one that avoids iterating over the potentially 
vast number of clusters not be evaluated, is to recognize 
that a valid cluster (containing a spin that is within Rc) 
of some heuristic strength may be built by adding a single 
spin to a valid sub-cluster of at least this same strength. 
This is formalized in the following theorem: 

Theorem 2 For any cluster, S, with a particular heuris- 
tic strength s and any spin in that cluster, i, there exists 
a sub- cluster C containing i of size ||C|| = ||<S|| — 1 with a 
heuristical strength c> s. 

This theorem follows from the heuristic strength defini- 
tion ensuring that S may be connected by interactions of 
strength s or greater and the following Lemma: 

Lemma 2 For any connected undirected graph, Q, and 
contained vertex, v, there exists a connected sub-graph, 
H, containing v but having one fewer vertex. 

Lemma 2 has a simple constructive proof. Starting from 
traverse Q until all vertices are visited; it doesn't mat- 
ter if any vertices or edges are traversed multiple times 
as long as we stop at the point at which all vertices have 
been visited. Since Q is connected, this must be possible. 
Because our stopping point comes once all vertices have 
been visited, the last vertex must have been visited only 
once. A sub-graph, with the desired property from 
Lemma 2 is formed by the vertices and edges of the tra- 
versed path excluding the final traversed edge and final 



visited vertex. This is illustrated in Fig. 25 

Given Theorem 2, we may find the clusters of size 
k with the highest heuristic strength by constructively 
building from clusters of size k— 1 with the highest heuris- 
tic strength. The algorithm works with three different 
lists: strongest clusters 5, potential clusters P, and the 
desired clusters D. The algorithm starts by adding into 
S all clusters of size one for each of the spins within the 
Rc cutoff and proceeds as follows: 

1. Take the strongest cluster C off of list P (which 
should be kept sorted) and add it to S. 

2. If D does not yet contain clusters of size k = 
||C||, add C to this D list as well. 
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FIG. 26: Schematic illustration of our algorithm to find 
the strongest clusters, the set of desired clusters D, with Nk 
quotas (in this example, 13 2-clusters, 7 3-clusters, and 2 
4-clusters). The dashed circle denotes the Rc cutoff. The 
black dot in the center represnent the central spin and other 
dots represent bath spins, (a) shows the initial population of 
1-clusters within Rc that must seed all selected clusters, (a)- 
(f) shows a progression of cluster selections, skipping a few 
selections at a time. The clusters grow from previously con- 
sidered clusters (which may are may not be in D depending 
on the Nk quotas) based upon the interactions strengths of 
the spin system being represented. 



3. Add into P any new cluster C (not already con- 
tained in S) that may be generated by extending C 
by one spin. 

4. Optionally, to minimize memory usage, remove all 
clusters that are too weak to be relevant (i.e., clus- 
ters that cannot compete or be built upon to com- 
pete for one of the Nk spots of strongest /c-clusters 
for any k). 

5. While \\D\\ < J2k^ki repeat from step 1. 

At the end of this process, D will contain, for each /c, the 
Nk strongest clusters of size k. A schematic example of 



this process is shown in Fig. 26 



Appendix B: Proof of Lemma 1 

We will prove Lemma 1 using induction and contra- 
diction. As the base case of our induction, note that La 



and 1/0 are constants, which follows from our prerequisite 
that L and L are constant when all coupling constants 
are taken to be zero. For convenience, let us denote the 
bath coupling power series as L s = fs({biesjes}) = fs 
and L s = fsdhesjes}) = fs- The b id parameters of 
fs and fs are occasionally dropped for convenience but 
are still implied (they should be regarded as power series 



functions of bij). By Eqs. ([5| and (J6|, 



fs = Yl 

CCS 

fs = fs/ n fe- 



rn 

(B2) 



CCS 



By induction, let us assume that fc({bi e cjec}) obeys 
the connected graph lemma for all C whose size is 
less than k. For the contradictory part of the proof, 
we assume that the lemma does not hold for some 
fs({biesjes}) with ||5|| = k. Thus, there exists in 
fs({biesjes}) some non-constant term, g({b ie s,jes}), 
and non-empty disjoint sets X,y C <S, such that 
g({bi e sjes}) does not depend on any bi e xjey- Such 
a term is therefore unaffected if we impose that all 
bizxjey De zero. By our assumed factorability property 
of L s = fs({bies,jes}), then 



fs({bies,jes})\ biexjey=0 = fxfy- 



(B3) 



By inductive reasoning, for all C C S such that CC\X ^ 

and c n y ^ 0, 



fc({hec,jec}) = const (B4) 



since all non-constant terms would contain bj e x,j ey fac- 
tors that are taken to be zero. From Eqs. (B2), (B3), 
and p4|, 



fs({h 



'iesjes 



}) 
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fxfy 
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(B5) 

Applying Eq. (Bl ) for S = X and 5 = 3^, the numerator 
and denominator above will cancel and we are left with a 
constant. This is a contradiction since this function, by 
the contradiction-proof assumption, should contain the 
non-constant g({bi e s,jes}) term that is not affected by 
imposing that b^xjey — 0- 
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